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Magneto-acoustic instability driven by cosmic rays and the waves that they 
excite was first studied by McKenzie and Webb and by Zank. In the present 
paper this instability is reconsidered. In general, the cosmic-ray plasma system 
is unstable. However, depending on how the time scale is chosen, the results and 
interpretation may be different from those of McKenzie and Webb and of Zank, 
particularly in the long- and short-wavelength limits. To illustrate this 
instability in an astrophysical context, we take our Galaxy and supernova 
remnants as examples. In those examples, length scales shorter than 2:5 kpe 
and 1:5 pe respectively are susceptible to this instability, and the growth times 
are of the order of 4 x 10° and 1:5 x 10* yr respectively. 


1. Introduction 


A common consensus about the interstellar medium (ISM) in our Galaxy is 
that the energy densities of various components, including thermal plasmas, 
turbulent plasmas, magnetic fields and cosmic rays, are roughly the same. Thus 
cosmic rays, which couple to thermal plasmas via embedded magnetic fields, 
play a dynamical role in various astrophysical objects. For example, study of 
cosmic-ray-modified shocks (see e.g. Drury & Volk 1981; Axford, Leer & 
McKenzie 1982; McKenzie & Volk 1982; Webb 1983; Duffy 1992) enables us to 
estimate the efficiency of shock acceleration of cosmic rays. A shock in a 
supernova remnant will accelerate cosmic rays, leading to modification of the 
shock, and hence affecting the evolution of the remnant (see e.g. Drury, 
Markiewicz & Volk 1989; Jones & Kang 1992). Modification of the solar wind 
termination shock by galactic cosmic rays or anomalous components has a 
direct effect on cosmic-ray modulation (see e.g. Ko, Jokipii & Webb 1988; Lee 
& Axford 1988; Ko 19915). Cosmic rays, regarded as an extremely hot tenuous 
gas, can support a galactic halo or drive a galactic wind (see e.g. Hartquist & 
Morfill 1986; Ko, Dougherty & McKenzie 1991; Ko 19910; Ipavich 1975; Zank 
19896; Breitschwerdt, McKenzie & Vólk 1991, 1993). All of these phenomena 
illustrate how cosmic rays are dynamically important in determining the 
structure of the ISM and the Galaxy. 

In this paper we discuss another aspect of the cosmic-ray-plasma system. We 
re-visit the work of McKenzie & Webb (1984) and Zank (1989a@) on cosmic-ray- 
driven magneto-acoustic instability. They considered a homogeneous or 
uniform system, which we may call a unidirectional wave system. In detail, the 
system comprises a thermal plasma in which a magnetic field is embedded, 
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cosmic rays and an Alfvén wave field propagating along the magnetic field. The 
wave field contains pre-existing and self-excited waves (those excited by 
cosmic-ray streaming instability). The plasma, cosmic rays and Alfvén waves 
are all regarded as fluids. A general conclusion is that the slow magneto-acoustic 
mode can be driven unstable when there is diffusion of cosmic rays. The energy 
feeding the instability comes from cosmic rays through the frictional force in 
generating Alfvén waves (for a detailed discussion see McKenzie & Webb 1984). 
In this work we consider a similar system to that in McKenzie & Webb (1984) 
and Zank (1989a), except that we take a wave-pressure-dependent cosmic-ray 
diffusion coefficient. Moreover, we allow for perturbation of the oppositely 
propagating wave. We confirm that in the linear regime the system is unstable 
in general, and show that only a single slow magneto-acoustic mode is unstable. 
(However, we find a new stable ‘island’, a set of measure zero in parameter 
space.) To further understand this unstable mode, we numerically investigate 
the relation between its growth rate and various parameters, which may give, 
say, the most unstable length scale. 

The hydrodynamic model for the cosmic-ray—plasma system with two Alfvén 
wave fields is presented in §2. A linear stability analysis is performed in §3, 
where we also analyse some asymptotic configurations, and show that the 
slower of the slow magneto-acoustic modes is unstable. Section 4 provides a 
detailed parametric study of the growth rate of the single unstable slow mode. 
A summary and a discussion of possible astrophysical manifestations of this 
instability is given in §5. 


2. The model 


The following hydrodynamic model describes, in a self-consistent manner, 
the interactions between a magnetized thermal plasma, cosmic rays and 
turbulent Alfvén waves. The Alfvén waves generated by cosmic-ray streaming 
instability act like a ‘clutch’ between the magnetized thermal plasma and 
cosmic rays. The model also includes the second-order Fermi effect. The 
governing equations are the ideal magnetohydrodynamic (MHD) equations 
augmented by energy equations for the cosmic rays and Alfvén waves (Ko 
1992; see also Drury & Volk 1981; McKenzie & Volk 1982). The overall mass 
continuity equation is simply that for the thermal plasma (since the cosmic rays 
and Alfvén waves are assumed to be massless fluids): 


0, 

"P V. (pU) =0, (1) 
ot 

where p and U are the density and flow velocity of the plasma. The overall 
momentum equation is 


dU 
pate. VU =—-V(E, E Fi PL +P) + Ix B+ pg. (2) 
Here P,, P, and P= are the pressures of the thermal plasma, cosmic rays, and 
forward- and backward-propagating waves respectively. B is the background 
magnetic field and J is the current density, and in the MHD approximation 
J = (Vx B)/u, (with u, the magnetic permeability of the vacuum). g is some 
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external acceleration (e.g. gravitation). The energy equations of the thermal 
plasma, cosmic rays, waves and background magnetic field are 


OE 
qt A =U. VE +P +P.) +I5.E+ pu.g, (3) 
OE, = v= P, 
E v.s = (vv, vn, (4) 
Okt y P, 
—wilVgt-U.VPtirx-— Vp. —— 
ot PUPUN otogr i VE, 27’ x 
OE s 
a tV.S=—J.E. (6) 
Here 
B 
am is (7) 
(Hop)? 


and the energy densities and fluxes are given by 


P 
HyeipU paper en 
g 


P 
E, = LE Ei-2Pi, E= 
and 
F, = U(jpU? + pe+P,), 


g 


F = (v.v. ue) - e. v E. 
v, +v 
(9) 
Fz =(ULV,) Ht +UP?, 


S=LExB, 
Ho 


where v, are the collisional frequencies of cosmic rays scattered by forward- and 
backward-propagating waves, and ¢ is the internal energy per unit mass of the 
plasma. In ideal MHD the electric field E = —U x B. The derivation of the 
cosmic-ray and wave energy equations (4) and (5) can be found in Ko (1992). 
(Note that the right-hand side of equation (28) of Ko (1992) lacks the term J. E, 
and a factor of 1 should be added to the right-hand side of equation (14) of that 
paper.) 

The set of equations is completed by Faraday's induction equation, which in 
the MHD approximation can be written as 


ôB 
z 7 YX (UB), (10) 
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along with the omnipresent constraint equation 
V.B=0. (11) 


It is clear from the energy equations (3)-(5) that energy exchange between 
thermal plasma, cosmic rays and waves is facilitated by two processes: the 
power of different pressure gradients (V ,. VP,, U. VP, and U.VP+), and the 
second-order Fermi effect (P./7). A detail discussion of this model can be found 
in MeKenzie & Webb (1984) (without second-order Fermi effect) and Ko (1992) 
(with second-order Fermi effect). 

To fix ideas, consider the cosmic rays as ultra-relativistic particles, with no 
cross-field-line diffusion. The diffusion tensor x and the growth time of the 
second-order Fermi effect 7 are then given by (Ko 1992) 

ce ViVa 


1 
FU DLE) prot79* 


v.v. VÀ 


v.v c Ca 
We take simple polytropic relations for the thermal plasma and cosmic rays (see 
(8)), in which y, and y, are constants. For ultra-relativistic cosmic rays y, = $. 
Moreover, we assume that the collisional frequencies are proportional to the 
Alfvén-wave pressures, i.e. v, = aPz (i.e. a represents the coupling between 
cosmic rays and plasma). 

Ignoring the spatial structure of the system (or considering a very localized 
spatial region), it can be shown that the second-order Fermi effect drives the 
system to a unidirectional wave system (Ko 1992). This essentially means that 
after a while only one of the waves survives. Now, let us turn on some spatial 
perturbation and ask what will happen. In the following sections we show that 
only in a set of measure zero in parameter space is the system stable or 
marginally stable. 


3. Linear stability analysis 


One of the first questions one would like to ask about a system like that 
described in 82 regards the nature of the simple equilibrium state and its 
stability. If the second-order Fermi effect or an external force is present, a 
simple uniform state (where state variables are independent of time and space) 
will not exist in the above system, (1)-(6), unless either forward- or backward- 
propagating wave vanish (see (12)) and g = 0. In a uniform state the absolute 
value of the plasma velocity is irrelevant, because the model is invariant under 
Galilean transformations. In this work we focus on uniform equilibrium states. 
Without loss of generality, we take the unperturbed plasma velocity and 
forward-propagating Alfvén-wave pressure to be zero (cf. McKenzie & Webb 
1984; Zank 19894). On the unperturbed uniform state yr, = (Po Us, Poo, Poo, P E0» 
B,) (with U, = 0 and Pj, = 0), we introduce a small perturbation dy = (dp, dU, 
ôP, ôP, ôP, ôB). After linearization, the relevant equations for the 
perturbation are as listed in appendix A with U, = 0, Pio = 0 and g = 0. To 
proceed, we make use of plane-wave perturbations (i.e. Fourier analysis) such 
as dy oc exp [i(k.x—wt)], where k is real and w may be complex. The system has 
two natural directions, B and k, and this leads to two distinct cases. In one B 
and k are parallel to each other while in the other they are not. However, after 
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some manipulation, the two cases give the same dispersion relation, which can 
be written as 


4V 4 )P. 
0- pn vi(n.- Wag t i vat) (v3 - K? Và 


x (V, + &V4)(V5 + LV a0 + iN? Vao) 


V?(1— p’) -> Vp the, 

2_y2 p n2 p2 lp 3/4" 40 

d |n: e y$-46 Vino P Vo A Vae 

= a? (Vp +3 HV 40)" | (13) 
re (Vp + V40)(Vn + HV 49 + inu? Vao) 


in which the phase speed of the plane-wave perturbation and the ‘sound’ 
speeds of thermal plasma, cosmic rays and waves are given by 
2 = Yao. P Yc Feo an = BP 00 


c0 ? wo 


Po Po 2p, 


respectively, p is the cosine of the angle between B and k (i.e. y = V 44. k/V4, k), 
and y is the normalized cosmic-ray diffusion coefficient, defined as 


(14) 


c? 


BaP os Vao (15) 


q= klp, lps 
We call 1, the diffusion length scale, and 1/y may be considered as the 
normalized coupling strength. 

As expected, the dispersion relation (13) includes that of ideal MHD. The 
terms in the second set of parentheses represent the Alfvén mode, and those in 
the square brackets (with a.o = aj) = 0, i.e. cosmic rays and waves are absent) 
represent the magneto-acoustic mode. The terms in the first set of parentheses 
represent the decaying mode of forward-propagating waves. (Of course, if we 
consider at the outset the unperturbed background to contain forward- 
propagating waves instead of backward-propagating ones, the decaying mode 
will then be that of the backward-propagating waves.) In the following we 
concentrate on the terms in the square brackets, and indeed we show that this 
compressive mode is unstable almost everywhere in parameter space. (This 
mode is shown to be compressive in appendix B. See also McKenzie & Webb 
(1984).) 

As one can imagine, it is not easy to thoroughly analyse the terms in the 
square brackets in (13), which involve six parameters (V,,, G59, à,,9, Geo, # and 
1). Nevertheless, one may proceed by studying some asymptotic configurations. 
We shall present four cases later. In any case, we know that the square 
bracketted expression has five roots (i.e. five modes) in general. McKenzie & 
Webb (1984) showed that the system is 'stable' in the short-wavelength limit 
(y >œ). They also showed that in the long-wavelength limit (7 = 0) the system 
is ‘unstable’ if a,, is sufficiently large. However, we must be very careful about 
the way in which they used terms ‘stable’ and ‘unstable’. What they (see also 
Zank 1989a) mean is that the imaginary part of the perturbation phase speed 
is zero in the short-wavelength limit, and is positive in the long-wavelength 
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limit together with sufficiently large plasma pressure. A mode is unstable 
(stable) if and only if Im (w) > 0(< 0). Therefore, even if Im (Vp) > 0, the 
system might be marginally stable (Im (w) = 0) in the long-wavelength limit 
(k = 0). On the other hand, even if Im (V,) = 0, the system might be unstable 
(Im (w) > 0) in the short-wavelength limit (4-00). 7 is a measure of the 
coupling between the plasma and cosmic rays. It depends on k and Pz. If we 
take k as constant (i.e. we consider a particular wavelength of perturbation), we 
find that the system is unconditionally stable when y = 0 (PZ —o0) and when 
3) -*o0 (Poo = 0). If we take Pgo as a constant (i.e. we consider a fixed diffusion 
length lp), we find that the system is unconditionally stable when 9 = 0 (i.e. 
k = 0), but unstable when 7 (i.e. k) is infinite (see $4 for details). Furthermore, 
we find that the system is stable (or marginally stable) under a combination of 
H, Vao aso and apo, When 

aso = 3(1 — 3H") V ho + 3(859)", (16) 
one of the five roots of the square bracketted expression in (13) is V, = — 14V4,, 
which says that this mode is marginally stable. The remaining four roots can be 
shown to have negative imaginary part (see appendix C), i.e. all modes are 
stable. Finally, there is an extra stable region, which is a consequence of the 
assumption of the absence of cross-field-line diffusion in our model. This region 
is ø = 0 (or the perturbation wave vector is perpendicular to the magnetic field). 
Indeed, the square of the phase speed is just the sum of the square of the 
‘sound’ speeds (V? = V4 tao +t (auo) --a2,). We believe that the system is 
stable only within these three (or four) regions of the parameter space, and they 
span a set of measure zero only. 

We now justify our assertion by presenting four asymptotic configurations: 
(i) small a,,; (ii) large 7; (iii) small 7; and (iv) a small deviation from the 
condition (16). Except that we are able to spell out which is the unstable mode, 
the result of the first three asymptotic analyses are basically the same as those 
of Zank (1989). For these first three cases we let 


V, = VO +eV 4- O(e), (17) 


where € = abo, 1/7 and 7 for cases (i), (ii) and (iii) respectively. Substituting (17) 
into the square brackets of (13), we find that VO satisfies ¥(V,) = 0 for cases 
(i) and (ii), and 2(V,) = 0 for case (iii), where 


Vpl — g? MAS 
diim Vio? rr Vi a ead pee bap (18) 
p p 
(V, ip VA)" 


G(V,) = F (Vp) tai (19) 


(V+ #V a0)" 
The first-order term is given by 


yo — Vp ipao [— Lip? Vgo/ VP + pVag)] (20) 
x Vy + UV go)? + (nu? 49)^] dF /dV ply 


for case (i), 


VO = ia? (21) 


ee a ps 
os Bg VA yo + LV 49) 


dV, 


(0) 
Vp 
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Ficure 1. Graphical solution of (18) and (19) (different zeroth-order dispersion relations). 
The solid, dashed and dotted lines correspond to the functions Z(V,)-- Vi —a2, Y(V,)+ 


30? 
V?—a?, and V? —a?, respectively. The figure clearly shows that Y(V,) may have three or five 
p— 0,0 p Pp M g y p 


real roots while F(V,) has four real roots. In this example a,, = a, = Vios Guo = 2V,,/5! and 
y = 08. 


for case (ii) and 


w = jq? 
VP = iaso 


p Val V) eet] d$ (22) 


(Vp) + VAS dV, vo 


for case (iii). In figure 1 the solid, dashed and dotted lines correspond to the 
functions F (Vp) + V$ —a$,, (Vp) + V5 — aso and V2, —a;, respectively. The inter- 
sections of the dotted and solid lines give the real roots of ¥(V,) = 0, while 
the intersections of the dotted and dashed lines give those of Y(V,) = 0. 
Consider case (iii). By moving the dotted line up and down in figure 1 (which 
corresponds to changing a2), we observe that 4 may have three or five real 
roots. In the case of five real roots, it is clear from the location of the roots and 
the derivative of Y at them that only one of the five V{ has positive imaginary 
part. In fact, this is the slower of the backward-propagating slow magneto- 
acoustic modes. In the case of three real roots the remaining two roots are 
complex conjugates (with real parts lying between 0 and —,4V,, i.e. the 
backward-propagating slow magneto-acoustie modes). Although one of these 
two modes has Im (V?) > 0, they are marginally stable because the zeroth- 
order growth rate vanishes when either k or P} is kept fixed (see the discussion 
above (16)). However, when we go to first order, the corresponding two Im (VQ) 
are complex conjugates. Im (V2?) is negative for the other three modes (the 
same argument applies as in the case of five real roots). Hence one and only 
one of the modes is unstable. Now consider cases (i) and (ii). Figure 1 shows 
that F has four real roots. Similarly to case (iii), the location of the roots 
and the derivative of F at them show that Im (V0?) > 0 for one of the modes 
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only. This is also the backward-propagating slow magneto-acoustic mode. 
We should point out that the fifth mode of case (i) or (ii) is stable. In fact, in 
case (i) this mode is V, = —4V4,— iyu? Vso d-O(a2,), while in case (ii) it is 
V, = — ip! Vas— HVao + O(3/1). 

The fourth asymptotic case stems from the condition (16). For fixed x and 
Vio, we consider a plane surface parallel to and at a small distance eV, from the 
surface defined by (16). As mentioned above, the system is stable when (16) is 
satisfied. Since we are only interested in the stability of the system in the 
neighbourhood of (16), it is sufficient to investigate the marginally stable mode, 


V, = —}4V 4, (the remaining four modes are stable, since they all have negative 
imaginary parts). After some manipulation, we obtain 
10215 pV. 754uV, a2,(—1+2inu) 
= —1uwV,)-e—— “+e 49 | A^ +2700 5 — — — 3 
V, UV 706—157 +e ul 2700 V4. x Apta) +O(e*), (23) 
where l 
M = {9V ho + aco +3 lano) + 8(1 — 4”) Vhoh/ Vico: | 


24 
N = {648[V3 +a + 3(a7)*] + 976(1 — 1) V4, — 6754? Và,/ V3, ) Mes 
Note that V,=—3V,, is the slower of the backward-propagating slow 
magneto-acoustic modes. (This can be demonstrated from figure 1, where the 
solid and dashed lines touch each other at V, = —$4V,.) 

These asymptotic results seem to suggest the system is unstable in general, 
and the unstable mode is the slower of the backward-propagating slow modes. 
However, when using a computer to explore the parameter space, we found that 
the slower mode is not necessarily the unstable one (see $4). Nevertheless, we 
are confident that the system is unstable almost everywhere, and there is only 
one unstable mode. In particular, it is one of the backward-propagating slow 
modes (of course, it will be one of the forward-propagating slow modes if we 
take Pao = 0 and Pio + 0 at the outset). Moreover, in most of the parameter 
space the unstable mode is the slower of the slow modes (see 84). 


4. Growth rates 


The asymptotic analysis in $3 indicates that the cosmic-ray-plasma system 
is stable only in exceptional cases. Furthermore, only one of the five magneto- 
acoustic modes is unstable. These features are confirmed by numerical 
exploration of the parameter space, taking account of the pressures of plasma, 
cosmic rays and Alfvén waves, the strength of the magnetic field, the angle 
between the magnetic field and the direction of propagation of the perturbation, 
the wavelength of the perturbation, and the coupling between cosmic rays and 
plasma. To facilitate discussion, we introduce several dimensionless quantities. 
We recall the projection factor u (— V ,,.k/V,, k) and the normalized diffusion 
coefficient y (= klp, which mingles the coupling «P, and the wavenumber k; 
see (15)). The various pressures are represented by the Alfvén Mach number of 
their ‘sound’ speeds (McKenzie & Webb 1984; Zank 19892): 


Oo Og (ano)? 
$ = s ; " = e " w = " 25 
is y oo uui ect E 
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Figure 2. The unstable mode is the slow magneto-acoustic mode ([Re (w/V,, k)| < #). The 
solid, dotted and dashed lines correspond to 7 = 0-01, 1 and 3 respectively. In the region 
above and to the right of each dividing line, the unstable mode is the faster of the two slow 
modes (i.e. the one with larger |Re («)]|), while in the region below and to the left of each line 
the unstable mode is the slower one. In this example f, = 0:1 and y = 31/2. 


Accordingly, we may take 1/V,, as the corresponding time scale. However, if 
we want to study the growth rates of different perturbation length scales, 1/v. 
(the collision time) or /5/V4, = K/V 3o = c*/3V*,v.. is more appropriate. We 
tentatively call /j/V4, the diffusion time scale. We define two dimensionless 
growth rates: 


| = Im (M,), PjzIm (72) = 7 Im (M,), (26) 


T.z 
k im V. 


w 
Vag k 
where M, is the Alfvén Mach number of the phase speed V, of the perturbation. 

The expression in square brackets in (13) is the dispersion relation of the 
magneto-acoustic mode of the cosmic-ray-plasma system. In general, the 
relation has five complex roots, and a positive imaginary part implies that the 
corresponding mode is unstable (except when 7 = 0, in which case the system 
may be marginally stable). We searched the parameter space quite extensively 
and found that, except for a set of measure zero (see $3), the system is unstable, 
and one and only one of the roots of the dispersion relation has positive 
imaginary part. Moreover, this unstable mode is one of the two backward- 
propagating slow magneto-acoustic modes (— 4 < Re (M,) < 0). (Of course, if 
we have Pz, = 0 while P7, + 0, the unstable mode is the forward-propagating 
one.) 

Each line in figure 2 separates the parameter space into two regions. In the 
upper/lower region the unstable mode is the faster/slower one (the one with 
larger/smaller [Re (M,)|) of the slow modes. It is also clear from figure 2 that as 
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Fievre 3. The growth rate T, (which is normalized to V,,k) as a function of various 
parameters. The solid, dotted and dashed lines correspond to fzo =0-1, 1 and 10 
respectively. (a) Bro = 01, Bao - 1, Bb E 1; (b) Bro = 1, Pao = 1, = 1; (c) Bro = 0*6, Boo = 1, 
y= 1; (d) =1, u= 1, fa = 1; (e) = 1, po 1, Bi. = 1. The gaps in (b)-(d) correspond to the 
condition (16). 
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FiGURE 4. (a) (P',)max, the maximum of I, with respect to 7, versus By. (b) (77)max» the value 
of 7 at which T, is a maximum, versus f. The solid, dotted and dashed lines correspond to 
MPa = 01, 1 and 10 respectively. In this example £. = 1 and w= 1. 


y increases, the region for the faster mode shrinks. Hence we confirm the claim 
that in most of the parameter space the unstable mode is the slower of the slow 
modes. 

We now present two sets of results on the growth rates. In the first we 
consider I',; thus V,,% is kept fixed and changing y corresponds to changing lp 
(or simply P;,). In the second we consider Ip; thus /5/V,, is kept fixed and 
changing 7 corresponds to changing k. (Note that V,, is always kept constant, 
since it is the normalization for other speeds.) Figures 3(a, b) show that when 
9 = 0 (or Po >), T, vanishes. In some cases (figure 3a), as 7 increases (or Pio 
decreases), I’, goes through à maximum and then tends to zero as 700 (or 

‘Py = 0). In other cases (figure 35), I, goes through a maximum first, then 
drops to zero and rises again to another maximum, and then tends to zero as 
9 00. We should point out that the diffusion coefficient is treated as a constant 
in McKenzie & Webb (1984) and Zank (1989a). They showed that for large 
enough P, the system is unstable (they mean Im (V,)-» 0) in the long- 
wavelength limit (7 = 0). In our case the diffusion length scale lp is inversely 
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Fiaure 5. Growth rate T, (which is normalized to V,,/1;) plotted as a function of various 
parameters. In this figure the solid, dotted and dashed lines correspond to f, = 0:1, 1 and 


10 respectively. (a) By = 1, Boo E lu = 1; (b) B» - 1, Boo = 1, n= 1; (c) 4 1, = 1, Ba = 
1; (d) n = 1, y = 1, £,) = 1. The gaps in (c) correspond to the condition (16). 


proportional to P;,, and the general trend of increasing P, is to stabilize the 
system. Hence, as y tends to zero, the system becomes less and less unstable, 
which is in contrast to their result. (Of course, we must remember that the 
model on which all these analyses based is invalid when the wave pressure 
becomes too large. Thus letting £,, tend to infinity is just for the sake of 
argument. Moreover, if we keep k and P}, fixed and let « tend to infinity, we 
recover the result of McKenzie & Webb (1984). Nevertheless, we should like to 
keep æa as a constant throughout this paper. The gaps in figures 3(b-d) 
correspond to the stability condition (16). In figure 4 (I,,) max is the maximum 
of T, with respect to y (or equivalently Pgo), and (5),,, is the 7 at which T, is 
a maximum. Note that there is a discontinuity in (7),4,, when £, satisfies the 
stability condition (16) (see figure 45). 

Figure 5 shows the behaviour of the growth rate Tp when the diffusion (or 
coupling) length scale lp (i.e. &P,,,) is kept fixed. Figure 5(a) illustrates once 
again that the growth rate is zero when 7 = 0. This simply comes from k = 0 
(the long-wavelength limit). Tp increases monotonically with y (i.e. k), and 
approaches a constant value (Tp) as 700 (the short-wavelength limit). This 
maximum value can be written down explicitly (see (21)): 


es ad 
E a prr acide || LY He 
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Figure 6. (a) (I',),., the value of I’, as >% and given by (27), versus £,. (I',),, is also 
the maximum of I’, with respect to y. (b) and (c) (9)g9 and (7)9;, the values of 7 when 
Dj/(Dj)J, — 09 and 0:5 respectively, versus 2, The solid, dotted and dashed lines 
correspond to fzo = 0-1, 1 and 10 respectively. In this example 2, = 1 and p = 1. 


where F is given by (18) and F (VO?) = 0. Note that McKenzie & Webb (1984) 
and Zank (1989a) concluded that the system is stable in the short-wavelength 
limit, because they discussed Im (w/V,, k) instead of Im (w) itself. Here we have 
Im (w) = Fp Vi, /lp. Keeping lp constant and letting 4 = klp tend to infinity 
gives a finite positive Im (v). Hence the system is unstable, indeed the most 
unstable, in the short-wavelength limit. The gaps in figure 5(c) come from the 
stability condition (16). Figure 5 shows that when the perturbation has length 
scales smaller than some critical scale, the system becomes very unstable. To fix 
ideas, we define, quite arbitrarily, two length scales (7), and (7), as the values 


2-2 
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of y when ['p/(I'p),. = 0:9 and 0:5 respectively. Figure 6 shows how (I'p),,, (1)o9 
and (5),, depend on f,,. These figures are of particular interest in astrophysics, 
since they give the major unstable length scales of the perturbation. We 
elaborate on this further in $5. 

Finally, numerical results show that the general trend of decreasing x is to 
reduce the growth rate (except for the small range around the stability 
condition (16); see e.g. figure 3c), i.e. to stabilize the system. 


5. Summary and discussion 

The instability of a system consisting of magnetized plasma, cosmic rays and 
self-excited Alfvén waves has been investigated using a hydrodynamic 
description of cosmic-ray propagation (Ko 1992). We have re-visited the 
cosmic-ray-driven magneto-acoustic instability discussed by McKenzie & Webb 
(1984) and Zank (1989a). We have considered a unidirectional wave background 
(e.g. Po + 0 and Pf, = 0). The second-order Fermi effect causes the perturbed 
oppositely directed wave (P$) to decay, but contributes nothing to the 
magneto-acoustic instability (recall the dispersion relation (13)). We have taken 
into account the dependence of the cosmic-ray diffusion coefficient on wave 
pressure, and hence deduced that the cosmic-ray-plasma system is uncon- 
ditionally stable in the long-wavelength limit (7 — 0), which is in contrast to 
the results of McKenzie & Webb (1984) and Zank (1989a) (see $4 and figures 
3a, b). However, this is not very important, since we knew from asymptotic 
analysis (see also Zank 1989a) and a numerical search that the system is 
unstable in almost all circumstances, and the stable region mentioned here is 
exceptional. The other stable regions are the short-wavelength limit (klp —oo 
with k fixed) and that where the condition (16) is satisfied. (Note that if we take 
the short-wavelength limit as kl, —oo with l, fixed then the system is unstable, 
in fact, the most unstable.) The source of free energy that fuels the instability 
comes from cosmic rays, and the energy is released through the frictional force 
in generating Alfvén waves when the cosmic rays stream through the 
magnetized plasma (for a detailed discussion see McKenzie & Webb 1984). 

A closer look at figure 4 reveals that when the length scale of the 
perturbation, k, is fixed, (T',) max (the maximum of the growth rate with respect 
to y when V,, k is fixed) attains its own maximum when f, S 0:5 and f,, « 1. 
A similar conclusion can be drawn when the growth rate is normalized to 
the collisional frequency, v_ = aP. Figure 6(a) shows that the largest of the 
local maxima of (L5), (the maximum of the growth rate for fixed aPpo 
occurs when 5-»00) occurs when f 1 and f, is very small. Indeed, 
analysis shows that for u = 1 and a given fpo (T5),, reaches its local maximum 


[Coolmax 
— _ Beo {[7 + 36 Bing + 48(5,,)*] fuo — [3 + 288,9 48(8 2, Puo t (5518 
8Bivo {[2 + 95, -12(,)*] — (7 + 128 20) [Broo + (Boa) T) 


at 1 
Bso = 1+3 Loo + 2Brco — 2Frool Boo + (Broo) T (29) 


This maximum is a monotonically decreasing function of fzo- When fj, tends 
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FiGURE 7. (a) (U,),, the value of Tẹ as y> (which is the maximum of P, with respect to 
7) versus zo- (b) and (c) (y), and (5y)y;, the values of y when T',/(T,),, = 0:9 and 05 


respectively, versus £7). The solid, dotted and dashed lines correspond to £, = 0:1, 1 and 10 
respectively. In this example, J. = 1 and p= 1. 


to zero, the maximum is proportional to (,,)? and f, tends to 1. With regard 
to this, we should make two remarks: (i) we cannot set f, equal to zero, 
because it is part of the normalization of the growth rate; and (ii) the absolute 
maximum of (I’p),, is located at £, = 0. To simultaneously explore the relation 
between the growth rate, the perturbation length scale and the wave pressure, 
we may introduce another dimensionless growth rate 


o 
* — Im (s. ico) wo^D 7 wo^k (3 ) 
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(However, the time scale c”/2ap, V^, is rather artificial.) Figure 7 (a) shows that 
for f,, = 1, (C4), has a maximum at f, + 02. (Note that (L«), — 8,9 as 
fus >.) After going through all these exercises, we conclude that when f, z 1 
and f,, is small, the system is very unstable. In other words, the cosmic- 
ray-plasma system is very unstable when some kind of equipartition occurs 
between the thermal plasma and magnetic field. Interesting enough, it is 
generally believed that such equipartition does happen in our Galaxy. 

The final outcome of this instability cannot, of course, be determined using 
the linear analysis presented in this paper. Numerical studies of the full system 
are therefore desirable. In any event, we can estimate the most unstable length 
scales, say, in our Galaxy. Following Ptuskin (1981), we take the large-scale 
magnetic field in our Galaxy as 3x 1075 G, and the gas number density as 
0:1 em^? (the warm intercloud region with temperature 10* K). The Alfvén 
speed is then 2 x 109 em s^!. The cosmic-ray diffusion coefficient in our Galaxy 
can be estimated by the abundances of secondary nuclei (e.g. Be and B) and the 
electron component (see e.g. the chapters Ptuskin and Dogiel in Berezinskii 
et al. 1990). A typical value is 10° cm? s~’. Therefore the diffusion length scale 
is about 10?* cm. Now, as mentioned above, the cosmic-ray-plasma system is 
very unstable when f, ~ 1 and f, is small (say 0:1-0-2). We can read off the 
most unstable length scales from figures 6 and 7. We concentrate on the more 
conservative estimate, (7)9.9. With the above parameters, figures 6(b) and 7(b) 
give (7)54 & 4. In other words, in our Galaxy the most unstable length scales of 
the cosmic-ray-driven magneto-acoustic instability are those smaller than A = 
2lp/ (1) & 2:5 kpc. The growth time is about 4 x 108/2, yr (figure 6a). We 
may take f.o © 1 in our Galaxy. McIvor (1977) showed that short-wavelength 
hydromagnetic waves can be damped efficiently by thermal conductivity, 
viscosity and ion-neutral collisions. In the above interstellar medium phase 
(i.e. the ‘hot intereloud regions’ in Melvor 1977), all damping rates are 
S 85 x 101*(27/A)* s"! (where A is the wavelength of the waves in cm). Thus 
for A z 2 pe dissipation by thermal conductivity, viscosity and ion-neutral 
collisions cannot suppress the instability described here. 

In the case of supernova remnants we take the magnetic field, density and 
diffusion coefficient as 107? G, 5 x 107? em? and 1075 cm? s^! respectively. The 
Alfvén speed is 3 x 10’ cm s^! and the diffusion length scale is about 3 x 10! cm. 
Therefore length scales S 1:5 pe are very unstable, and the growth time is about 
1:5 x 10!/f, yr. 

Finally, we should like to end this paper with two questions. The first is about 
the stability condition (16). What is so special about the combination of 
pressures such that the system becomes (marginally) stable, at least in the 
linear regime? A partial answer may be provided by the fact the component 
ôP, of the eigenvector of the marginally stable mode vanishes (see (B 2) in 
appendix B). A natural conclusion is that dP, = 0 is a necessary condition for 
stability. Of course, this does not resolve the puzzle of that specific combination. 
If the perturbation becomes so large that the linear analysis presented here is 
no longer valid, will the condition (16) still hold ? The second question concerns 
the evolution of a system containing both forward- and backward-propagating 
waves. If the system is spatially homogeneous, the second-order Fermi effect 
drives it towards a unidirectional state (Ko 1992). However, the analysis here 
shows that this state is unstable to spatial perturbation. So what is the final 
state? 
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Appendix A 


In this appendix we give the linearized equations for the general system (cf. 
McKenzie & Webb 1984): 


9» Lu, ps pV. 8U 4-8U. V p, 4-6p V. U, — 0, (A 1) 
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2R LU, VàB— B. V9U--B, V.9U-- 9U. VB, —9B. VU, - BV. U, = 0, 
. . . < $ (A 6) 
along with the linearized constraint equation 
V.ôB = 
and ôV 4 is given by = E (A 7) 
V P V. (A 8) 


a (Ho Po)? 2Po 


Appendix B 

Here we show that the mode represented by the terms in the square-brackets 
in the dispersion relation (13) is compressive, i.e. V2 = dPiota,/dp. First we write 
down some of the components of the eigenvector corresponding to the 
eigenvalues of the square-bracketted expression in (13): 


ôP, = a, dp, (B 1) 
Vo +347, 
ôP, = a3, — e 40 — Sy, B2 
eT Y ERAS E Vaa id 
z - a Vo tiu, Vaol Vp +34V a0) | 
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Vo- 8B = - (mp FEV 49-8 e (A VÀ, óp, (B 4) 
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Recall the linear perturbation of magnetic pressure: ee 
1 Po i 
OP, = —B.6B = |= ] V,,.6B. (B 6) 
Ho 0. 


Making use of the square bracketted expression in (13), we obtain the linear 
perturbation of the total pressure of the system 


OP, = OP, + OP, + OP, + OP, = V? dp. (B 7) 
Thus we have shown that the mode is compressive. Since ôP}, does not appear 
in the square brackets in (13), and in fact it belongs to the decaying mode (the 


first set of parentheses in (13)); hence it is independent of p. Thus we ignore dP}, 
in the above calculation. 


„Appendix C 


In this appendix we consider a special case to illustrate that the imaginary 
part of V, is less than or equal to zero when the condition (16) is satisfied. If (16) 


MHD instability driven by cosmic rays 41 


holds then the square-bracketted expression in the dispersion relation (13) 
becomes 


V, -iuV, 2 
0 = EET [0 ea a Me 
V, + iuV, V,—&uV, 
—a?, e #40 _ 41 —y?) Vi, 2 * 49), (C1) 
d Vo UVa tinue? Vao eae V,— LV 40 l 
Take the simple case p = 1. The above equation is further reduced to 
Vo +3Va0 1 fn \2 2 V, +3V a0 
een —$V49) ^3 — 0n Te |. 2 
0 VF Vae (V, + Vao) Vp —3V49) —$(2,) — aco V, Vo t iyVa (C 2) 
The root V, = —1V,, is real, and in the following we show that the remaining 


three roots have negative imaginary parts. First let V, = X -- Y, with real X and 
Y. The square-bracketted expression in (C2) becomes two simultaneous 
equations 


X? + 3Vigg X° — [BY? + 29V gq Y + $(859)? + aco] X 
BY? +3 9V a0 Y HiV ho Hildu) +3460] Vso = 9, (3) 
(3Y + 9V4q) X* + (3Y +39 V a0) Vao X 
{Y° 4 Vao Y? + Glano) + aco] Y + NVa Vao Hilaw) h= 0. (C 4) 
Eliminating X from these equations, we obtain a horrendously complication 
equation for Y: 


0 = (3Y + V4) (27 648 F? + 829445 V,, Y° 
+ 3456(12a2, + 16(0,,,)? --9V*, +287? V?) Y? 
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+ 24V,,,[16(3a2, + 4(a%5)?)(21a2, + 34(a,,)?) 
+324V2, ,(5a2, + 8(0,,)?) + 729V4, 
+ 65? V2,(180a2, + 544(a,,,)? + 243V*.,) + 72% V*,] Y* 
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t 39? V4,[16(261a75, + 100802 (an0)? + 640a,5) 
+ 216(37a2, + 60(a,,,)?) V4, 4- 3645V*,] 
+ 10874 V*,(32a2, + 192(a7,.)? +81V%,,)} Y? 
+ 9Vao{4[16(3a5, + 4(a,,)*)? 
+ 108(9a5, + 24475 (2,,,)" + 4(,,,)*) Và + 729029 Vir) 
t 39? V2,,[16(45a4, + 33602, (aj9)* + 128(a,,,)*) 
t 216(11a2, + 12(a,,,)?) V4, + 729V 40] 
+ 367° V44(664(a2,,)* + 21V4,)) Y* 
305, 9° Vial (12055 + 16 (a7)? + 9V2,)(12a7, + 64 (azo)? + 27 V0) 
+ 121? V4,(64(2,,)* + 27V 40) Y 
+ 72085 9° Vio (8(8,,)* + 3V5)]. (C 5) 
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Although this equation looks terribly difficult to deal with, we found that all the 
coefficients are positive. It is clear that Y must be negative in order to satisfy 
(C 5). Hence the imaginary parts of all four roots of (C 2) are less than or equal 
to zero. Thus the system is stable when the condition (16) holds. The one for 
general p is orders of magnitude more complicated than (C 5) (it has terms up 
to Y?? instead of Y?), but the same conclusion can be reached. 
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